Thermodynamics of the in-medium chiral condensate 
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The temperature dependence of the chiral condensate in isospin-symmetric nuclear matter at 
varying baryon densities is investigated starting from a realistic free energy density of the correlated 
nuclear many-body system. The framework is thermal in-medium chiral effective field theory which 
permits a systematic calculation of the quark mass dependence of the free energy density. One- and 
two-pion exchange processes, virtual A(1232)-isobar excitations and Pauli blocking corrections are 
treated up to and including three-loop order. It is found that nuclear matter remains in the Nambu- 
Goldstone phase with spontaneously broken chiral symmetry at least in the range of temperatures 
T < 100 MeV and baryon densities up to about twice the density of normal nuclear matter. 
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The chiral condensate (qq), i.e. the expectation value 
of the scalar quark density, plays a fundamental role as 
an order parameter of spontaneously broken chiral sym- 
metry in the hadronic low-energy phase of QCD. The 
variation of (qq) with temperature and baryon density 
is a key issue for locating the chiral transition boundary 
in the QCD phase diagram. The melting of the conden- 
sate at high temperatures and/or densities determines 
the crossover from the Nambu-Goldstone phase to the 
Wigner-Weyl realization of chiral symmetry in QCD. 

It is thus of principal interest to perform a systemati- 
cally organized calculation of the thermodynamics of the 
chiral condensate. Such a calculation requires knowledge 
of the free energy density as a function of the quark mass 
(or equivalcntly, as a function of the pion mass). The ap- 
propriate framework for such a task is in-medium chiral 
effective field theory with its explicit access to one- and 
two-pion exchange dynamics and the resulting two- and 
three-body correlations in the presence of a dense nuclear 
medium. 

Previous studies of the in-medium variation of the chi- 
ral condensate were mostly concerned with the density 
dependence of (qq) at zero temperature, using differ- 
ent approaches such as QCD sum rules lj or models 
[H [3] based on the boson exchange phenomenology of 
nuclear forces. Temperature effects have been included 
in schematic Nambu - Jona-Lasinio (NJL) approaches 
[H [5] . Such NJL models work with quarks as quasipar- 
ticles and provide useful insights into dynamical mecha- 
nisms behind spontaneous chiral symmetry breaking and 
restoration, but they do not properly account for nucle- 
ons and their many-body correlations, a prerequisite for 
a more realistic treatment. 

The present work extends a recent chiral effective 
field theory calculation [5] of the density-dependent in- 
medium (qq) condensate to finite temperatures T. A 
related chiral effective field theory approach to nuclear 
matter at T = performing rcsummations to all orders 
has been reported in ref.[7j- Corrections to the linear 



density approximation are obtained by differentiating the 
interaction parts of the free energy density of isospin- 
symmetric nuclear matter with respect to the (squared) 
pion mass. Effects from one-pion exchange (with m r - 
dependent vertex corrections) , iterated l7r-exchange, and 
irreducible 27r-exchange including intermediate A(1232)- 
isobar excitations, with Pauli-blocking corrections up to 
three-loop order are systematically treated. The domi- 
nant nuclear matter effects on the dropping condensate 
are supplemented by a further small reduction due to 
interacting thermal pions. To anticipate the result: we 
find that the delayed tendency towards chiral symmetry 
restoration with increasing baryon density p, observed at 
T = [H] in the same framework, gets gradually softened 
with increasing temperature. An approximately linear 
decrease of the quark condensate with increasing p is re- 
covered at temperatures around T ~ 100 MeV. However, 
no rapid drive towards a first order chiral phase transi- 
tion is seen, at least up to p < 2 p where po = 0.16 fm -3 
is the density of normal nuclear matter. 

Consider the free energy density T = pF(p,T) of 
isospin-symmetric (spin-saturated) nuclear matter, with 
F(p, T) the free energy per particle. In the approach to 
nuclear matter based on in-medium chiral perturbation 
theory [5J [5] the free energy density is given by a sum of 
convolution integrals of the form, 
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where K,i,K,2 and /C3 are one-body, two-body and three- 
body kernels, respectively. The last term, the so-called 
anomalous contribution A(p, T) is a special feature at 
finite temperatures |10] with no counterpart in the cal- 
culation of the groundstate energy density at T = 0. 
As shown in ref. [5] the anomalous contribution arising in 
the present context from second-order pion exchange has 
actually very little influence on the equation of state of 
nuclear matter at moderate temperatures T < 50 MeV. 
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denotes the density of nucleon states in momentum space. 
It is the product of the temperature dependent Permi- 
Dirac distribution and a kinematical prefactor p/2ir 2 
which has been included in d(jp) for convenience. Mm 
stands for the (free) nucleon mass. The particle density 
p is calculated as 





FIG. 1. Three-loop contributions to the free energy density 
of nuclear matter that are relevant for the in-medium chiral 
condensate. Left diagram: pion self-energy correction; cen- 
tral and right diagram: two-pion exchange Hartree and Fock 
terms involving the nnNN contact interaction proportional 
to the low-energy constant c\m\. 
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This relation determines the dependence of the effective 
one-body chemical potential p(p, T; Mm) on the thermo- 
dynamical variables (p, T) and indirectly also on the nu- 
cleon mass Mm- The one-body kernel K\ in eq.(l) pro- 
vides the contribution of the non-interacting nucleon gas 
to the free energy density and it reads [5] 
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The first term in K\ gives the leading contribution (den- 
sity p times nucleon rest mass Mm) to the free energy 
density. The remaining terms account for (relativistically 
improved) kinetic energy corrections. 

Our starting point is the Feynman-Hellmann theo- 
rem which relates the temperature dependent in-medium 
quark condensate (qq) (p, T) to the quark mass deriva- 
tive of the free energy density of isospin-symmetric (spin- 
saturated) nuclear matter. Using the Gell-Mann-Oakes- 
Renner relation m 2 f 2 = — m g (0|gg|0) one finds for the 
ratio of the in-medium to vacuum quark condensate 
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where the derivative with respect to m 2 . is to be taken 
at fixed p and T. The quantities (0|g<7|0) (vacuum quark 
condensate) and f„ (pion decay constant) are to be un- 
derstood as taken in the chiral limit, m q — > 0. Likewise, 
m 2 stands for the leading linear term in the quark mass 
expansion of the squared pion mass. 

In the one-body kernel K\ the quark (or pion) mass 
dependence is implicit via its dependence on the nucleon 
mass Mm- The condition dp/dM^ — applied to eq.(3) 
leads to the following dependence of the effective one- 
body chemical potential p. on the nucleon mass Mm- 
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The nucleon sigma term um = (N\m q qq\N) = 
m 2 8Mm /dm 2 measures the variation of the nucleon 
mass Mm with the quark (or pion) mass. Combining 



both relationships leads to the m^-derivative of the one- 
body kernel: 
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In the limit of zero temperature, T = 0, the terms in 
eq.(7) reproduce the linear decrease of the chiral conden- 
sate with density. The kinetic energy corrections account 
for the (small) difference between the scalar and the vec- 
tor (i.e. baryon number) density. In the actual calcu- 
lation we use the chiral expansion of the nucleon sigma 
term <jm to order 0(m 4 ) as given in eq.(19) of ref.[5]. The 
empirical value of the nucleon sigma term (at the physi- 
cal pion mass = 135 MeV) is a N = (45 ±8) MeV [TT] , 
In numerical calculations we choose the central value. 

The two- and three-body kernels, K-2 an d £3, re- 
lated to one-pion exchange and iterated one-pion ex- 
change have already been given in explicit form in ref. [8] . 
Their derivatives with respect to the squared pion mass, 
dK-2^/dm 2 , are hence obvious and do not need to be 
written out here. The same applies to the anomalous 
contribution A(p, T) (see eqs.(14,15) in ref.[8]) and to 
the two- and three-body kernels related to 27r-exchange 
with excitation of virtual A(1232)-isobars (see section 6 
in ref. [9]). In case of the one-pion exchange contribu- 
tion we include the m^-dependent vertex correction fac- 
tor r(m Tr ) as discussed in section 2.1 of ref. [5]. The short- 
distance contact term which produces a T-independent 
correction of order p 2 to the in-medium condensate is 
treated exactly in the same way as in ref. [6], i.e. all terms 
with a non-analytical quark-mass dependence generated 
by pion- loops are taken into account. 

Let us now turn to some three-loop contributions which 
are new and of special relevance for the in-medium chiral 
condensate. The first one comes from the pion selfenergy 
diagram shown in Fig. 1. It gives rise to the following 
derivative of the two-body kernel: 
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the 7T7T low-energy constant £3 ~ 3. 
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The chiral ttttNN contact vertex proportional to cim 2 
generates 27r-exchange Hartree and Fock diagrams, also 
shown in Fig.l. Concerning the free energy density 
pF(p, T) or the equation of state of nuclear matter their 
contributions are actually almost negligible. However, 
when taking the derivative with respect to m 2 as re- 
quired for the calculation of the in-medium condensate, 
these contributions turn out to be of similar importance 
as other interaction terms. The corresponding contribu- 
tion to the derivative of the two-body kernel reads: 
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with the auxiliary function: 
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The 27r-exchange Hartree diagram with one cimj-vertex 
contributes the following piece to the three-body kernel: 
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while the three-body term associated with the 2ir- 
exchange Fock diagram with one ci to 2 -vertex gives: 
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Last not least we incorporate the effects of thermal 
pions. Through its m 2 -derivative the pressure (or free 
energy density) of thermal pions gives rise to a further 
reduction of the T-dependent in-medium condensate. In 
the two-loop approximation of chiral perturbation theory 
including effects from the 7T7r-interaction one finds the 
following shift of the condensate ratio in the presence of 
the pionic heat bath [T2Tll4j : 
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with the functions Hi t3 (m n /T) defined by integrals over 
the Bosc distribution of thermal pions: 
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We proceed with a presentation of results. As input we 
consistently use the same parameters in the chiral limit 
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FIG. 2. Density dependence of the chiral condensate in 
isospin-symmetric nuclear matter taken at a temperature 
T — 100 MeV. Starting from the linear density dependence 
(dashed curve) characteristic of the free nucleon Fermi gas, 
the following interaction contributions are successively added: 
one-pion exchange Fock terms (l7r), second order (iterated) 
pion exchange, irreducible two-pion exchange (2tv), two- and 
three-body contributions from 27r exchange with intermediate 
A(1232)-excitations (A), and two-pion exchange with tttvNN 
vertex involving the chiral low-energy constant cimj. Pauli 
blocking effects are included throughout. 



as in our previous works [6], namely: f n = 86.5 MeV, 
g A = 1.224, ci = -0.93 GeV- 1 and M N = 882 MeV. 
Concerning the contact term representing unresolved 
short-distance dynamics, we recall from ref. [6j that its 
quark mass dependence, estimated from recent lattice 
QCD results for the nucleon-nucleon potential [TS], is 
negligibly small compared to that of the intermediate and 
long range (pion-exchange driven) pieces. 

It is worth pointing out that in-medium chiral pertur- 
bation theory with this input produces a realistic nuclear 
matter equation of state [9], including a proper descrip- 
tion of the thermodynamics of the (first-order) liquid-gas 
phase transition. Apart from temperature T, the addi- 
tional "small" parameter in this approach is the nucleon 
Fermi momentum pp in comparison with the chiral scale, 
^frr ~ 1 GeV. Our three-loop calculation of the free en- 
ergy density is reliable up to about twice the density of 
normal nuclear matter. It can be trusted over a tem- 
perature range (up to T ~ 100 MeV) in which the hot 
and dense matter still remains well inside the phase of 
spontaneously broken chiral symmetry. 

Fig. [2] shows a representative example, at T — 100 
MeV, displaying stepwise the effects of interaction contri- 
butions to the density dependence of (qq) (p, T) from the 
chiral two- and three-body kernels £2,3. As in the T = 
case studied previously [6], the pion-mass dependence of 
correlations involving virtual A(1232)-excitations turns 
out to be specifically important in delaying the tendency 
towards chiral symmetry restoration as the density in- 
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FIG. 3. Ratio of chiral condensate relative to its vacuum value 
as function of baryon density p in symmetric nuclear matter, 
for different temperatures up to T = 100 MeV as indicated. 
The reduction effect due to thermal pions is included. 

creases. Once all one- and two-pion exchange processes 
contributing to K-2 and /C3 arc added up, the chiral con- 
densate at T = 100 MeV recovers the linear density de- 
pendence characteristic of a free Fermi gas. However, 
this recovery is the result of a subtle balance between 
attractive and repulsive correlations and their detailed 
pion-mass dependences. Had we taken into account only 
iterated one-pion and irreducible two-pion exchanges, the 
system would have become instable not far above nor- 
mal nuclear matter density as seen in Fig. [2] In fact, 
this instability would have appeared at even much lower 
densities in the chiral limit (m ff — > 0). This emphasizes 
once more not only the importance of terms involving 
A(1232)-excitations, but also the significance of explicit 
chiral symmetry breaking by small but non-zero quark 
masses in QCD and the resulting physical pion mass, in 



governing nuclear scales. 

Fig. [3] shows the systematics in the variation of the chi- 
ral condensate with temperature T and baryon density 
p. These results include all nuclear correlation effects 
and also the (small) additional shift arising from thermal 
pions. The latter correction is visible only at the high- 
est temperature considered here (T = 100 MeV) where 
the chiral condensate at zero density begins to deviate 
from its vacuum value. According to recent QCD lattice 
simulations the actual crossover transition at which the 
quark condensate (qq) drops continuously to zero occurs 
around a temperature of T c ~ 170 MeV [16] . 

At zero temperature, the hindrance of the dropping 
condensate at densities beyond normal nuclear matter 
comes primarily from three-body correlations through /C3 
which grow rapidly and faster than K2 as the density 
increases. The heating of the system reduces the influ- 
ence of K.3 relative to JC 2 continuously as the temperature 
rises, so that their balance at T = 100 MeV produces a 
small net effect in comparison with the free Fermi gas. 

In summary, this is the first calculation of the quark 
condensate at finite temperature and density that sys- 
tematically incorporates chiral two-pion exchange inter- 
actions in the nuclear medium. Correlations involving 
intermediate A(1232)-excitations (i.e. the strong spin- 
isospin polarizability of the nucleon) together with Pauli 
blocking effects are demonstrated to play a crucial role 
in stabilizing the condensate at densities beyond that of 
equilibrated nuclear matter. The results reported here 
set important nuclear physics constraints for the QCD 
equation of state at baryon densities and temperatures 
that are of interest e.g. in relativistic heavy-ion colli- 
sions. In particular, we find no indication of a first order 
chiral phase transition at temperatures T < 100 MeV and 
baryon densities up to about twice the density of normal 
nuclear matter. 
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